Evolutionary Biology

The general random walk model (GRW) of Hunt (2006) is used to infer directional evolution in mean trait values from sparse fossil data by modeling phenotypic change as the accumulated result of small steps with mean step sizes and step variances. Using simulations and real data cases, Ergon (2026) showed that the step variances can be estimated reasonably well only when the mean trait values have small measurement errors, while for fossil data with realistic measurement errors they appear to be extremely difficult to find, and they are often found to be negative. In the simulations Ergon (2026) assumed that the true phenotypic mean values were known. Here, I essentially repeat these simulations under the assumption that only mean trait values with large measurement errors are known, and based on weighted mean squared error (WMSE) comparisons the conclusion is that weighted least squares (WLS) is a better method than GRW. A second conclusion is that WLS is a better method also in the possibly rare cases with large measurement errors where the GRW parameters are estimated well. The GRW method is simply not flexible enough to handle such cases. A third conclusion is that Akaike Information Criterion (AIC) results for GRW models with large measurement errors relative to the step variance may be overly optimistic.

Sharks exhibit extraordinary morphological diversity across a wide range of ecological niches, yet large-scale, high-resolution digital datasets of their internal anatomy remain limited. Here we present an open-access 3D shark anatomical repository derived from published X-ray computed tomography (CT) data, featuring manually segmented and systematically annotated models of the chondrocranium, visceral arches, axial skeleton, musculature, and viscera in standard STL format. The dataset comprises 117 individuals, representing 72 species across 25 families and all nine extant shark orders, with 115 full-body reconstructions and two head-only models. This open-access dataset offers a comprehensive resource for comparative anatomy, biomechanical simulations, evolutionary developmental biology and biomimetics research of extant sharks.

The origin of novel phenotypes can influence access to new ecological resources, which may have positive, neutral, or negative effects on subsequent phenotypic diversification. In this study, we tested the macroevolutionary consequences of a pair of putative functional innovations occurring in deep-sea fishes of the order Stomiiformes. Integrating phylogenetic comparative methods, micro-CT scans, and external body measurements, we recover a mosaic of diversification trends associated with these innovations. We found some evidence for elevated evolutionary rates in tooth morphology associated with the predatory dragonfishes, which possess a gap between their vertebral column and skull that exposes the notochord and enables neck-like flexibility. However, a second novelty building upon the first, a functional neck joint enabling extreme cranial kinesis, was linked to faster rates of skull evolution. Our results suggest that innovations that help shift ecological roles and overcome functional constraints related to those roles, like gape-limitation in prey depauperate habitats, may play an important role in promoting phenotypic diversification. This work builds on a growing body of evidence highlighting how the deep sea promotes phenotypic diversity, generating the extreme forms that are celebrated by scientists and the public alike.

Brain size and structure are hypothesized to be adaptively designed to satisfy the behavioral requirements of securing food and living socially. The importance of these socioecological and sociobiological selective forces in brain evolution is constantly debated. Socioecological divergence is striking in the Neotropical ant genus Neoponera: N. apicalis is a generalist solitary predator forming small colonies of ~100 whereas N. commutata colonies are approximately 10 times as large and workers pheromonally organize cooperatively raids only on Syntermes termite colonies. We interspecifically compared the size and structure of the compound eyes, size and number of antennal glomeruli, mosaic brain scaling and synaptic processing (microglomeruli-MG). Our results indicate that N. apicalis workers have a larger number of ommatidia, antennal lobe glomeruli, and allometrically larger antennal and optic lobes than N. commutata. These sensory traits were associated with differences in higher-order processing architectures in the mushroom body (MB) microglomeruli (MG). N. commutata workers had an allometrically larger MB, perhaps due to their socially complex chemical foraging communication, although MG density in N. apicalis was higher in both the MB lip and collar, regions associated with processing olfactory and visual information, respectively. The increase in MG density in N. apicalis may be associated with higher demands for navigation, learning, and memory, as well as a higher density of antennal lobe glomeruli to support prey odor discrimination. In contrast, N. commutata workers had larger ommatidia and antennal lobe glomeruli. Larger ommatidia correlate with their diurnal/nocturnal habits and a larger MB Our findings indicate that differences in behavioral performance demands associated with socioecological differentiation are reflected in variation in visual and olfactory system structure, brain size, mosaicism, and synaptic organization. Our results support both social and ecological brain hypothesis as drivers of mosaic brain evolution.

Untargeted metabolomics offers a powerful lens for quantifying high-dimensional phenotypic variation within and among species in nature, but has yet to be widely adopted in evolutionary ecology. Some important initial questions are whether metabolome composition differs among populations, and to what extent such variation is genetic or plastic. Here, we use untargeted liquid chromatography tandem mass spectrometry to characterize the relative abundance of 5,939 molecular features of the threespine stickleback (Gasterosteus aculeatus) liver metabolome. Native lake populations differ in metabolome composition, reflecting effects of sex, size, geography, and population ecotype (benthic versus limnetic). Stickleback from these lakes were translocated to found new populations in nine recently fishless lakes, permuting fish ecotypes across benthic and limnetic lake habitats. Several generations later, metabolomes in these experimental populations reflect effects both of their genetic ancestry (e.g., taurocholic acid, a cholane steroid bile acid, was elevated in limnetic-ancestries), as well as their present habitat (e.g., acylcarnitines). Additionally, ecotypes transplanted into a habitat to which they were maladapted exhibited a distinctive metabolomic profile. We conclude that stickleback exhibit both heritable and plastic among-population differences in liver metabolome, which could represent an important phenotypic basis of rapid evolution, population divergence, and perhaps local adaptation.

Color polymorphism can facilitate local adaptation, maintain intraspecific diversity, or reflect early stages of speciation. Clownfishes (Amphiprion spp.) typically display a simple black, orange, and white pattern, but the saddleback clownfish (Amphiprion polymnus) shows striking variation in melanism and the number of vertical bars, which are thought to play a role in species recognition. In 2024, a revision on iNaturalist split A. polymnus into multiple species based solely on color pattern and geographic range. This raises the question of whether these morphs represent true species or intraspecific polymorphism, which we tested using genomic and image-based data. We sampled 97 individuals from seven populations across the species range and quantified color patterns from standardized photographs. Phenotypic and genomic analyses reveal a complex pattern of divergence. Image analysis identified three distinct phenotypic clusters, with A. polymnus, A. annamensis, and A. laticlavius each showing consistent differences in saddle shape and vertical bar extent. ADMIXTURE resolved three distinct genetic groups corresponding to the morphs. Pairwise FST (0.54-0.71) and dxy indicate extremely high differentiation between A. polymnus and A. annamensis, consistent with species-level divergence, whereas A. laticlavius shows much lower differentiation from A. polymnus (FST 0.09-0.18) and higher differentiation from A. annamensis (FST 0.64-0.66). Overall, phenotypic and genomic data show structured variation, but the status of A. laticlavius remains ambiguous. Our study reveals clear and structured divergence across the full range, yet the taxonomic interpretation of this variation remains inherently challenging. The key question remains: do these patterns reflect a single polymorphic species or a complex of closely related species?

The preparation of multi panel figures remains a labor intensive step in scientific publication. Albeit there are specific tools available to solve this problem, they are often highly specialized, difficult to install, or time consuming to learn. Griphus is a standalone graphical application designed for rapid composition and experimentation with multi panel figures, developed by and for zoological taxonomists. Functions specifically designed for multi panel composition include automatic figure numbering and placement, aspect ratio operations, spacers, layout rotation, layout suggestions, and automatic generation of figure legends, including scale bar descriptions. The software can perform both spatial interpretation of images on the canvas and work with a simple, editable layout formula. It also enables instant multi panel composition, with numbered images and automatic contrast selection for the numbers, obtained simply by loading images. User defined parameters such as target printable dimensions, resolution, spacing, and color mode are preserved throughout the work. The program produces coordinated outputs consisting of the final composite figure, a readable file describing the layout structure, and a .gri file storing images, transformations, and parameters for exact regeneration. Griphus is intended as a complementary tool to professional image software, providing a simple and efficient environment for constructing high quality multi panel figures.

Humans and other vertebrates contain two estrogen receptors (ERs), ER-alpha and ER-beta, which mediate the physiological actions of three estrogens: estrone (E1), estradiol (E2) and estriol (E3). Of these three estrogens, in vivo, E2 is the strongest transcriptional activator of ER-alpha and ER-beta, E1 is next most active, followed by E3. We studied transcriptional activation of human ER-alpha and ER-beta by E2, E1 and E3 in African green monkey kidney (COS-7) cells, which we compared with studies of estrogen stimulation of ER transcription in human em-bryonic kidney (HEK-293) cells. To our surprise, in COS-7 cells, E3 had the lowest half-maximal response (EC50) for human ER-alpha and ER-beta than either E2, which was second most active estrogen, or E1. In contrast, for human ER-alpha and ER-beta transfected into HEK-293 cells, E2 was the most active estrogen, followed by E1 and E3. Similar results were found in COS-7 cells and HEK-293 cells transfected with elephant shark ER-alpha and ER-beta. Thus, under some conditions, E3 is a more active estrogen than either E2 or E1. This suggests that E3 may be a novel physiological ligand for the ER in some mammalian cells.

Coral reefs are constructed by colonial cnidarians whose survival depends on the coordinated growth and physiological integration of thousands of interconnected polyps. While coral skeletons have been extensively studied, the internal three-dimensional organization of coral tissues remains poorly resolved, limiting our understanding of how reef-building corals function as integrated modular organisms. In this study, we established a contrast-enhanced X-ray tomography (XRT) workflow for decalcified coral tissues, enabling detailed visualization and quantitative comparison of internal polyp architecture across four reef-building species with distinct colony forms: Acropora cervicornis, Acropora millepora, Montipora capitata, and Pocillopora damicornis. Importantly, this methodology resolved previously inaccessible patterns of tissue organization and structural connectivity among neighboring polyps. The two Acropora species shared a conserved axial - radial organization but differed in mesenterial morphology, whereas M. capitata exhibited complex, entangled mesenterial networks that connected both neighboring and distant polyps. In contrast, P. damicornis displayed superficial connectivity restricted to the coenosarc. Together, these results suggest that internal tissue architecture is an evolutionarily flexible trait, shaped by ecological and developmental pressures rather than strictly by shared ancestry. Our XRT workflow thus provides a new comparative framework for understanding how corals function as integrated living colonies.

Eusocial Hymenoptera present the highest known recombination rates among metazoans, which evolved several times independently among bees, ants and wasps. Several hypotheses have been proposed to explain this observation, including stronger selection for recombination caused by coevolving parasites and pathogens, and strong sexual selection among haploid males due to male-biased sex ratios among reproductive individuals. In this article, we explore the effects of haplodiploidy and differential selection between sexes on the evolution of recombination, by analyzing a three-locus model in which selection for recombination stems from negative epistasis between selected loci. Our analytical predictions are compared with the results of individual-based simulations in which deleterious mutations occur along a linear chromosome. Our results show that, at mutation-selection balance for deleterious alleles, increasing the strength of selection against deleterious alleles (due to the effect of male haploidy and/or sexual selection) tends to reduce selection for recombination. However, an increase in the overall magnitude of negative epistasis (which may also be due to male haploidy and/or sexual selection) combined with the fact that recombination only occurs in females may increase selection for recombination substantially. Our model also shows that, in conditions favoring recombination, increasing recombination in meioses leading to parthenogenetic ovules (and male offspring) may yield stronger benefits than in meioses leading to fertilized ovules (and female offspring).

Allele frequency (AF) timeseries allow us to directly observe the dynamics of evolution at a genetic level. However, extracting useful inferences from AF timeseries has proved difficult due to the model uncertainties and noisiness inherent in AF change at fine temporal scales. Here we present three new permutation tests --- which do not assume a model of evolutionary change or a parametric statistical model --- to detect AF timeseries features of evolutionary interest. The features identified by these approaches are: 1) any evolutionary change (as opposed to apparent change due to measurement error); 2) directional selection; 3) fluctuating selection with a propensity to change sign (negative autocorrelation). We are not aware of existing tests for features 1 and 3. Feature 2 is commonly tested using standard evolutionary models such as the Wright-Fisher; we show that the permutation approach has comparable statistical power. We apply our new approaches to AF timeseries data from D. melanogaster and D. pulex.

Observed group sizes rarely match the size that would maximize what each member gets from belonging. We propose a two-part theory in which group size is regulated by two related conflicts: insider-outsider conflict over admission, and within-group conflict as crowding, competition, and social tensions intensify with size. Three strategies are available: admission, exclusion, and fission. The first part shows that even when exclusion is unavailable, fission dynamics alone drive group size away from the optimum in both directions, with the pattern set by how prospective joiners encounter groups and by the geometry of fission. When joiners compare groups across a shared landscape and fission is asymmetric, the standing distribution is bimodal: supra-optimal large groups coexisting with a sub-optimal mode of small groups, the pattern characteristic of fission-fusion societies. The second part promotes exclusion and fission to active decisions: incumbents weigh the per-capita cost of accommodating entry ({beta}) against the costs of coordinated exclusion (c +{gamma} N*) and fissioning (F). A single inequality, {beta} > c +{gamma} N*, partitions populations into two regimes: where it holds, exclusion is viable and groups lock at the optimum size; where it fails, groups grow past the optimum and cycle through recurrent fission. Modal group size, fission frequency, and exclusion behavior together identify which regime governs a population -- a set of predictions applicable across fishes, social insects, birds, and mammals including primates and human foragers.

A recent study (Manjarrez et al., 2026) showed that the classification of cortical dendritic spines into stubby, thin, and mushroom subtypes is unstable under rotation. That result criticizes the categorical scheme but leaves an open question. What is the actual structure of spine morphology once the viewing angle is controlled? Here we answer it. We analyzed 228 spines from layer II pyramidal neurons in the H01 nanometer-resolution reconstruction of human temporal cortex. We first quantified the source of instability. We found that rotating dendritic segments by 90 degrees about their axes shifted the apparent spine height and head width in opposite directions across the population, thereby confirming orientation-dependent measurement error. Furthermore, to obtain measurements free of this artifact, we developed the Spine Morphometry Hub (SMH), a 12-point anatomical landmark framework that characterizes each spine in all three orthogonal planes and extracts geometric, voxel-based, and mesh-based metrics. All morphometric distributions were unimodal and right-skewed. Density-based clustering assigned most spines to noise, and a Monte-Carlo test against a discrete two-type null model confirmed that this pattern is incompatible with categorical subtypes. We also confirmed that apical and basal spines were statistically indistinguishable. Unlike previous reports of a spine continuum, all based on orientation-dependent measurements, our framework removes the viewing-angle confound itself, so the continuum we observe cannot be attributed to a projection artifact. Hence, our framework will be useful to quantify dendritic-spine remodeling in neurological disorders, in which spine shape has long been observed but never measured against an orientation-invariant morphometric standard. HighlightsO_LISpine Morphometry Hub (SMH) measures spines free of viewing-angle error C_LIO_LISMH was validated as an orientation-invariant morphometry framework C_LIO_LIRotating dendrites by 90{degrees} shifts spine height and head width oppositely C_LIO_LIAll morphometric distributions are unimodal and right-skewed, not categorical C_LIO_LISMH could be used to quantify dendritic-spine remodeling in neurological disorders C_LI

Mammal tails have long been recognized for their diversity of morphological form and function, however, there remains a substantial gap between the motivation to understand and emulate the various performance functions of the tail and what is known about tail anatomy. In this study, we were motivated to discover the anatomical foundations of the fast, whipping motions of the tail of the lesser Egyptian jerboa (Jaculus jaculus), which may aid in the quick changes of direction as the animal escapes from predators using ricochetal bipedal hopping. We employed microCT scans, dissections, and museum data to describe the musculoskeletal anatomy of the jerboa in comparison with the laboratory mouse (Mus musculus) and rat (Rattus norvegicus). While many aspects of tail anatomy are conserved across these species, the jerboa does possess unique characteristics such as an extremely long tail arising from caudal vertebral elongation, development of extensive dorsal musculature differentiated into lateral and medial components to increase points of skeletal attachment, and a novel anatomical feature - the bi-lobed cranial transverse process - which serves as a supernumerary dorsal tendon attachment site and possible brace to protect the ventral tendons and intrinsic muscles for a section of caudal vertebrae which likely experiences high mechanical stress.

Many proposed mechanisms for the evolution of cooperation among unrelated individuals rely on relatively demanding cognitive abilities that are not widespread across taxa. In contrast, individual heterogeneity is a pervasive feature of animal groups, encompassing differences in personality as well as physical and cognitive traits. Such heterogeneity can promote the evolution of cooperation, yet its role has received comparatively little attention, particularly as a source of variation giving rise to social organization such as leadership. A specific form of leadership can emerge under unstable environmental conditions, when some individuals become better suited than others to initiate action and influence the behavior of their peers. Unlike fixed dominance hierarchies, emergent leadership can rapidly adjust to changing environmental conditions, thereby reshaping group organization. Because it does not require the maintenance of stable hierarchies, this form of leadership can arise even in species that do not have the cognitive capabilities to sustain complex social structures. In this work, we investigate the combined effects of individual heterogeneity and emergent leadership on the evolution of cooperation using an evolutionary game-theoretic model in which individuals may assume the roles of leaders or followers according to their strength, representing individual differences in suitability to prevailing environmental conditions. We examine different levels of population heterogeneity together with increasingly complex strategy sets requiring progressively greater informational requirements, allowing individuals to condition cooperation on their own strength, leadership role, or both. Our results show that the interplay between leadership and heterogeneity promotes the evolution of cooperation, particularly when only a small fraction of individuals act as leaders. Under these circumstances, cooperation evolves even when individuals employ the simplest possible strategies. Under harsher ecological conditions, cooperation can be sustained by more sophisticated strategies, specifically by conditional strategies that prescribe cooperation when individuals are strong or leading and defect when acting independently. Author summaryIn this study, we propose that emergent leadership mediated by individual diversity can boost the evolution of cooperation in animal groups. Building on growing evidence on the heterogeneity of animal capabilities and personalities, we focus on the fleeting leadership that emerges in animal groups when facing rapidly changing environmental conditions. We suggest that this type of leadership that emerges from individual differences in strength--a generic quality encompassing those characteristics that make an individual more fit to lead in a given situation--does not require complex cognitive capabilities from the animals and represents a valid alternative to more demanding strategies proposed in the past to explain the evolution of cooperation. Using an evolutionary game theory model, we show that if a population includes a few strong players, these can become influential leaders and guide the actions of their peers to achieve cooperation. Although the naive strategy of always cooperating is sufficient for cooperation to evolve, the introduction of more complex strategies leads players to cooperate only when they are more likely to be recognized as influential leaders. These strategies are more effective in promoting cooperation under unfavorable ecological conditions and are also more robust against exploitation by defectors.

The number of chromosome sets per nucleus is a fundamental trait, but why this number is nearly always two for multicellular eukaryotes remains unclear. Chromosomes are made of nucleic acids, which possess abundant phosphorus (P). Therefore, producing new chromosomes, as well as generating new cells and organismal growth, demands substantial phosphorus. Yet, because P is often limiting in nature, P availability could influence the prevalence of diploidy versus polyploidy. Here, we compare growth rates of diploid and triploid Potamopyrgus antipodarum, a freshwater snail, relative to P availability. Because diploid P. antipodarum are obligately sexual while obligately asexual individuals are polyploid, costs associated with sensitivity to P limitation in polyploids could also help explain the maintenance of sexual P. antipodarum. We raised juvenile diploid and triploid snails on either P-adequate or P-deficient diets and found that independent of P availability, juvenile triploid asexual snails grew faster and harbored higher P content as adults than sexual diploid conspecifics. Together, these results suggest life-history advantages of polyploidy or asexual reproduction that exacerbate rather than ameliorate the cost of sex. These outcomes suggest that P availability is unlikely to be a main driver of ploidy polymorphism or the maintenance of sex in P. antipodarum.

Genomic offset (GO) is increasingly used to predict population maladaptation risk under climate change, with larger offsets assumed to indicate greater vulnerability. Despite rapid adoption in conservation planning, it remains unclear how sensitive GO estimates are to key methodological choices, including SNP set composition, genotype-environment association (GEA) methods, and the specific GO metric used. Empirical validation against observed population dynamics also remains limited. Here, we evaluate the methodological robustness and predictive performance of GO using multidecadal demographic monitoring data from Lycaeides butterflies, a system with short generation times and high fecundity that may facilitate rapid adaptive responses. GO estimates were broadly consistent across SNP sets, regardless of composition or size, with climate-associated and randomly selected SNPs yielding largely concordant values. Consistency across GEA methods was moderate and depended on the SNP set used. In contrast, GO metrics differed substantially in the magnitude of maladaptation estimated, suggesting they capture distinct biological signals and should not be treated as interchangeable. Crucially, GO was a poor predictor of observed population trends, regardless of SNP set composition, GO metric, or GEA method, both at sites used to fit GEA models and when extrapolated to independent demographic sites. These findings suggest that, while GO provides a valuable conceptual framework for assessing potential maladaptation, its quantitative estimates and predictive power are sensitive to methodological choices and species-specific biological context. We therefore urge careful alignment of GO metric assumptions with conservation objectives, along with rigorous empirical validation, before GO estimates are used to inform management decisions.

Empirical studies have demonstrated that mutation rates may change with individual condition, such as in the case of stress-induced mutagenesis. This has led to the hypothesis that condition-dependent (or "plastic") mutation rates could be selectively favoured, as the increased production of new mutants in times of maladaptation enhances evolvability, the ability to undergo adaptive evolution. However, while empirical evidence for condition-dependent mutation rates is accumulating, theoretical models studying their evolution are lacking. Here, we employ an individual-based simulation approach to examine the evolution of condition-dependent mutation rates in a changing environment. We find that condition-dependent mutation rates consistently evolve when the environment changes at an intermediate pace. Furthermore, populations with condition-dependent mutation rates are substantially better adapted to their (changing) environment. Finally, the evolutionary dynamics of condition-dependent mutation rates are both accelerated and destabilised when the mutation rate is self-referential (i.e., when mutator loci affect their own mutation rate). We conclude that condition-dependent mutation rates (and thus evolvability) can readily evolve in changing environments. Significance statementMutation provides the raw material for evolution. Mutation rates thus tune evolvability, the ability to undergo adaptive evolution: if mutation rates are too low, evolution is impeded; if mutation rates are too high, adaptive traits cannot be maintained. Using a theoretical model, we explore the evolution of plastic mutation rates that systematically depend on the condition of the organism and its environment. An example is stress-induced mutagenesis in bacteria, which is implicated in the evolution of antibiotic resistance. We show that plastic mutation rates readily evolve, providing "well-timed" variation specifically when organisms are poorly adapted. Such plastic mutation rates thus facilitate better adaptation to changing environments, and their evolution provides an example of evolvability itself evolving.

Recent evidence of extensive call sequence use in non-human primates has led to the theory that syntax evolved to mitigate the constraints of their genetically fixed repertoires, before vocal production learning later emerged in humans. However, evidence of similarly extensive sequence repertoires in an open-ended vocal production learner--the Western Australian magpie (Gymnorhina tibicen dorsalis)--offers a unique opportunity to explore potential alternative pathways to syntactic communication. Our previous work revealed fledgling magpies learn group-specific repertoires of structured call sequences from their social contacts, with more sociable individuals acquiring larger repertoires earlier in development. Notably however, the individual vocal segments that combine to form their calls and call sequences were shared across groups and emerged as early as the first week post-fledging--suggesting the underlying vocal elements may not be learned. Here we utilised acoustic neighbourhood-based dimensionality reduction to compare clustering patterns of vocal segments across magpie fledgling developmental stages, and between fledglings and adults. We found no evidence of acoustic development over time, and no significant distinction between fledgling and adult productions of the same vocal segments. The same coarticulatory effects--where a vocal element is produced differently when combined with another--and geographic variation established previously in adults were supported in fledglings too. These findings support that the vocal building blocks underpinning magpie call sequences are innate, suggesting usage learning better explains how fledglings learn to combine calls. In a species capable of open-ended production learning, this suggests learning to combine existing signals may be more adaptive than productively learning new ones. Rather than evolving solely to compensate for genetically fixed repertoires, syntax may have evolved as a flexible, convergent solution to the various challenges of expanding communicative capacity--whether due to genetic constraints, cognitive limitations or the cost of establishing new meaning in novel signals.

The global study of freshwater clams in the genus Corbicula is frequently confounded by invasive androgenetic lineages that experience mitochondrial DNA capture and clonal propagation. In contrast, the endemic Corbicula of Sulawesi's ancient lakes reproduce sexually, offering a uniquely reliable system for mitochondrial population genetics. This study provides the first population-genetic framework for two endemic species, Corbicula possoensis (Lake Poso) and C. linduensis (Lake Lindu), using the cytochrome c oxidase subunit I (COI) marker. We analysed 90 newly generated COI sequences from C. possoensis (six stations) and C. linduensis (three stations), integrated with reference sequences from GenBank, to assess genetic diversity, population structure, and phylogeographic patterns. Hierarchical AMOVA revealed deep divergence between the two lakes ({Phi}_CT = 0.607), consistent with prolonged independent isolation rather than a single shared vicariance event, as the two species do not form a sister pair in the phylogeny. Within Lake Poso, C. possoensis exhibited exceptionally high genetic diversity (24 haplotypes; h = 0.876; {pi} = 0.016) and pronounced micro-geographic structuring into three phylogeographic zones (North: Tentena and Siuri; East: Tando Nceppo and Busogo Beach; Southwest: Bancea and Pendolo), each characterised by distinct haplogroups. Remarkably, the maximum divergence between zones (K2P = 2.33%) approached the interspecific distance between C. possoensis and C. linduensis (K2P = 2.42%), indicating that within-lake mitochondrial divergence has reached near-interspecific levels. Conversely, C. linduensis displayed near-panmixia and extreme genetic depauperation (3 haplotypes; h = 0.246; {pi} = 0.0004), indicating long-term demographic stasis within a restricted habitat. The deep phylogeographic zonation in C. possoensis suggests that its discrete populations should be treated as separate Management Units (MUs) in conservation planning to preserve locally adapted gene complexes, whereas the severely depauperate gene pool of C. linduensis renders it critically vulnerable to environmental disturbance and invasive species, warranting urgent IUCN Red List assessment. To validate these mitochondrial boundaries and inform future conservation strategies, multi-marker and genome-wide reassessments are strongly recommended.